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Abstract 

Decaying electron magnetohydrodynamic (EMHD) turbulence in three dimensions is studied via 
high-resolution numerical simulations. The resulting energy spectra asymptotically approach a A;~^ 
law with increasing Rb, the ratio of the nonlinear to linear timescales in the governing equation, 
consistent with theoretical predictions. No evidence is found of a dissipative cutoff, consistent 
with non-local spectral energy transfer and recent studies of 2D EMHD turbulence. Dissipative 
cutoffs found in previous studies are explained as artificial effects of hyperdiffusivity. In another 
similarity to 2D EMHD turbulence, relatively stationary structures are found to develop in time, 
rather than the variability found in ordinary or MHD turbulence. Further, cascades of energy 
in 3D EMHD turbulence are found to be suppressed in all directions under the influence of a 
uniform background field. Energy transfer is further reduced in the direction parallel to the field, 
displaying scale dependent anisotropy. Finally, the governing equation is found to yield a weak 
inverse cascade, at least partially transferring magnetic energy from small to large scales. 

PACS numbers: 
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I. INTRODUCTION 



Turbulence plays a crucial role in a wide variety of geophysical and astrophysical fluid 
flows. In this paper we present results on a speciflc variety of plasma turbulence in which 
the flow consists entirely of electrons, moving through a static background of ions. The 
equation governing the electrons' self-induced magnetic field is then 

^ = - V X [J X B] + R-^' V^B, (1) 

where J = V x B, and Rb = crBo/nec, with a the conductivity, Bq a measure of the field 
strength, n the electron number density, e the electron charge, and c the speed of light. 
See for example lo|, who derived this equation in the context of magnetic fields in the 
crusts of neutron stars. More generally though, it is applicable in many weakly collisional, 
strongly magnetic plasmas, so other applications could include the Sun's corona or the 
Earth's magnetosphere. 

Turbulence governed by ([1]) is known as electron MHD (EMHD), Hall MHD, or whistler 
turbulence. Based on its (at least superficial) similarity to the vorticity equation governing 
ordinary, nonmagnetic turbulence, 

^ = V X [u X w] + Re'^ VV, (2) 



where now w = V x u, 



10| argued that ([T]) would initiate a turbulent cascade to small 



lengthscales, thereby accelerating neutron stars' magnetic field decay beyond what ohmic 
decay acting on large lengthscales could achieve. They suggested in particular that the 
turbulent spectrum would scale as k^^, with a dissipative cutoff occurring at A; ~ Rb- 
However, recent work has highlighted the fundamental differences between equations ([1]) 



and 



21 



Wareing and HoUerbach 2009], henceforth referred to as WH09). In ([2]) the 
dissipative term contains more derivatives than the nonlinear term, so on sufficiently short 
lengthscales the dissipative term will always dominate and hence there is a dissipative cutoff. 
In contrast, in ([T]) the two terms both contain two derivatives, so it is conceivable that the 
nonlinear term will always dominate, even on arbitrarily short lengthscales. [21„ WH09] 
found precisely this effect in decaying 2D EMHD turbulence. 

In this paper we present high-resolution numerical simulations of ([1]) in a three- 
dimensional (3D) periodic box geometry, designed specifically to address such questions 
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as to whether there is a dissipative cutoff or not, and whether the couphng is local or not. 
In contrast to previous 3D simulations , Biskamp et al. 1999, Biskamp and Miiller 

1999, Cho and Lazarian 2004]), we do not employ hyperdiffusivity, which disrupts this fea- 
ture that the two terms in ([T]) have the same number of derivatives, and hence introduces 
an artificial dissipative cutoff ([21, WH09]). These previous simulations found no difference 
between 2D and 3D turbulent spectra, with both having a /c"^/^ scaling. In our recent study 
of the governing equation in two dimensions with normal diffusivity ([21, WII09]), we found 
no evidence for a dissipative cutoff, and the turbulent spectrum scaled as A;~^/^, broadly 
consistent with previous results. We now consider whether 3D EMHD turbulence displays 
the same scaling characteristics. Finally, we also consider the question of whether ([T]) is 
capable of yielding an inverse cascade in 3D (we use hyperdiffusivity for this set of runs). 



II. NUMERICS 



We solve ([I]) by treating the z- independent parts of B as before ([2l|, WH09]). For the 
z-dependent parts we expand and By in triple Fourier series in x, y and z, with B^ 
then given by V ■ B = 0. Time integration is achieved through a second order Runge- 
Kutta method, with the diffusive term treated exactly. We employ standard pseudospectral 
techniques for the evaluation of the nonlinear terms, with dealiasing according to the 2/3 
rule. The code employs the MPI library 22] and the FFTW library ([8, Frigo and Johnson 
2005]) to achieve massive parallelisation on a suitable supercomputer. We performed a 
variety of runs, typically employing 64 processors, with the highest extending to /c = 170 in 
Fourier space, corresponding to = 512 collocation points in real configuration space. Due 
to the two derivatives in the nonlinear term, the required timesteps are unfortunately very 
small, roughly proportional to 1/(A^^). Values as small as ~ 1 x 10^^ were used, requiring 
0(2 X 10^) timesteps to reach t = 0.2. 



A. Initial Conditions 

Since our interest is in freely decaying, rather than forced, turbulence, we need to carefully 
consider the nature of our chosen initial conditions. We will present results for three different 
sets of runs. 
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First, to study homogeneous forward cascades, we start off with random 0(1) energies 
in all Fourier modes up to A; = {k^ + fc^)^''^ = 5, making sure that the energy is evenly 
distributed between the three components. After initialisation the overall amplitude of the 
field is rescaled to ensure that the rms value of |B| = 1 at t = 0. 

Second, to study nonhomogeneous forward cascades, we start off with the same initial- 
isation as above, but now add a uniform field Csx, where C = 1, 2, 4 or 8. This field is 
simply added to the component directly, resulting in a suitably modified equation ([1]). 

Third, to explore the possibility of inverse cascades, we return to the C = case without 
any large scale magnetic field and now inject energy into modes in the range 10 < k < 20. 
The question then is how much of this initial energy moves to k < 10, and how much moves 
tok> 20. 

Finally, for all three sets of results, each individual run was repeated with a number 
of different random initial conditions, to ensure that the results presented here are indeed 
representative. 



B. Ideal Invariants 



Equation ([T]) has some useful associated diagnostics, corresponding to quantities that are 
conserved in the ideal, R~b^— > 0, limit. Specifically, in three dimensions we have equations 
for the energy, 

f;-jB'.W^-Ri'j3^-.W, (3) 

and the magnetic helicity. 



d 1 
~dt2 



A ■ B dy = -R^^ / B • J dV^, (4) 



where A is the vector potential, defined by B = V x A. Note though that in the presence 



Berger 1997]), let alone 



of a uniform background field, helicity is not even defined 
conserved. 

Whereas in 2 dimensions, we also had the additional quantity of the mean squared mag- 
netic potential or anastrophy, this is not conserved in 3 dimensions. This quantity is partic- 
ularly important for the inverse cascade in 2D EMHD, where such a cascade is thought to 
be driven by a forward cascade of energy and an inverse cascade of anastrophy ([18, Shaikh 
and Zank 2005]). It is unclear then, whether an inverse cascade will occur in 3D. 



In addition to the physical insight that these various integrated quantities yield into 
the nature of the Hall nonlinearity, they also offer useful diagnostic checks of the code. 
Reassuringly, we found that both of them (except helicity in a uniform field) were satisfied 
to within 0.25% or better by all of our runs. 



III. RESULTS 



A. Large-scale initial conditions 

A characteristic statistical quantity of a turbulent system is the energy spectrum. 3D 
hydrodynamic turbulence exhibits the spectrum Ek oc k^^^^, the famous Kolmogorov law 



l4 . Kolmogorov 1941]). In MHD turbulence, the energy transfer is altered by the Alfven 



effect 
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15l . Iroshnikov 1964, Kraichnan 1964]), leading to a flatter energy spectrum 



Ek (X k Recent studies of 2D EMHD turbulence have found, via methods which all 



employ hyperdiffusivity, a 5/3 Kolmogorov spectrum for small scales /ccL > 1 



equivalent to 



18l . Biskamp 



k > 0{Rb), and a steeper 7/3 spectrum for longer wavelengths 

et al. 1996, Biskamp et al. 1999, Dastgeer et al. 2000, Dastgeer and Zank 2003, Cho 

and Lazarian 2004, Shaikh and Zank 2005]). Of the studies which considered 3D EMHD 
fl 

turbulence, [3[ numerically found a scaling of 7/3 for kde < 1 consistent with a local spectral 
energy transfer independent of the linear wave properties. However, the authors of that 
paper noted the scaling in 3D was only marginally verified, since the short extent of the 
inertial spectral range was dominated by the bottleneck effect, causing local enhancement of 
the spectrum above the inertial range power law at the point of transition from the inertial 
range to the dissipation range. The more abrupt the transition (or the higher the degree of 
hyperdiffusivity), the more pronounced the effect. 

In the left plot of Figure [U, we show the energy spectra of our solutions for Rb = 
10, 30, 100 & 300, evolved to a time t = 0.2. The energy spectra have been stationary 
since approximately t = 0.16 and time averaging between 0.16 and 0.2 reveals an identical 
spectrum and no further information. We interpret this to mean our simulations are resolved 
and evolved to a suitable time for inspection of the now quasi-stationary cascade. It is worth 
noting at this point that 3D EMHD turbulence takes a longer time to reach quasi-stationarity 
than 2D EMHD turbulence, which achieves such a phase by t ~ 0.1. 
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Fourier power spectra 
1 • • • ^ 




Compensated power spectra 




FIG. 1: Energy spectra of homogeneous 3D EMHD turbulence at t = 0.2. On the left are shown 
spectra for Rb = 10, 30, 100, & 300. On the right, compensated spectra k'^Ek for the same range 
of Rb. 

The energy spectra all start out much the same at low k, peak around k = 3 and then 
lower Rb spectra smoothly drop off with increasing k whilst higher Rb spectra maintain 
a linear gradient in the log-log plot. Transfer of energy to higher k is then more efficient 
at higher Rb, with Rb = 300 having a scaling of ~ A;~^/^. The spectra are asymptotically 
approaching an energy spectrum oc k~'^, where we propose ~ 2. In the right plot 
of Figure [Tj we show compensated energy spectra to show this approach to k'^Ek = 1 
with increasing Rb- Our value of u is the same as that predicted by lo!] for 3D EMHD 
turbulence. Their prediction was calculated using a phenomenology based on Kraichnan's 
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also considered a theoretical 



arguments (the whistler effect) and is reproduced in 
prediction neglecting the whistler effect and derived a scaling of 7/3. Simulations have shown 
the whistler effect has little effect on the energy spectrum of 2D EMHD turbulence, but it 
would seem that it may have a considerable effect on 3D EMHD turbulence. 

Regarding the bottleneck effect, we see no evidence for it in the spectra. This is not a 
great surprise as we see no sign of a dissipative cutoff, which causes the effect. By definition, 
the dissipation scale should occur when the local value of Rb is 0(1) in equation ([1]). It is 
unclear though when this occurs since the definition of Rb does not involve length scales. 
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If the coupling is purely local in wavenumber, then this definition does involve length scales 
after all, since the Bq that should be used is the field at that wavenumber only, rather than 
the total field. That is, according to the definition of ll| where this argument was first 
developed, we have 

R's = Rb{B'/B) (5) 

where the primed quantities are the small-scale local values and the unprimed the large-scale 
global. If we now suppose a energy spectrum, then B' / B ~ and so R'^ is reduced 
to 0(1) when k ~ Rb- So, at Rb = 100 we would see a dissipative cutoff at A; ~ 100, which 
we do not. It is possible to reconcile the situation by realising that this argument crucially 
depends on the coupling being local in Fourier space: if this does not hold then R'^ = Rb 
and there is simply no definite dissipation scale. 

In agreement then with our 2D study of EMHD turbulence, the nonlinear term is able 
to dominate at all length scales and the coupling is therefore non-local in Fourier space. 
Again, hyperdiffusivity has previously masked the effect of the nonlinear term at high k. It 
has also introduced the bottleneck effect which has disguised the true scaling of the energy 
spectra. With normal diffusivity, we believe we have performed the first true simulations of 
3D EMHD turbulence. 

In Fourier space then, EMHD turbulence continues to bear a strong resemblance to 
ordinary MHD turbulence. We would like to know if this resemblance carries over into real 
configuration space. In Figure [2l we show slices through the three component field datacubes 
at t = 0.1. Across the top row, we show a slice of the B^ field at x = 0.5 on the left, a 
slice at y = 0.5 in the middle and a slice at z = 0.5 on the right. Across the middle row 
we show the same slices for By and across the bottom row for B^. Large numbers of small, 
independent vortices dominate the fields, characteristic of fully developed turbulence. 

In Figure [3l we show the same fields at t = 0.2. It is clear that the fields resemble those 
at t = 0.1, implying 3D EMHD turbulence is much more structured than classical and MHD 
turbulence, where a fully developed turbulent field would bear no resemblance to the initial 
field. This appears to be a unique characteristic of decaying EMHD turbulence, in both two 
and three dimensions. 

We would also like to address the energy decay of the field, with particular respect to any 
dependency of the decay rate on the value oi Rb- ^ reported that the energy dissipation 
rate is independent of the value of the dissipation coefficient, represented hj Rb here. In 
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FIG. 2: Plots of the Rb = 300 solution in real configuration space at t = 0.1. For full details of 
the slices through the datacube, see the text. 



contrast, 
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WH09] found the energy decay is strongly dependent on the Rb parameter, 
with the decay rate proportional to Rb^- We find the same behaviour here, as shown in 
Figure IH 



B. Large-scale initial conditions in the presence of a background field 



3D EMHD turbulence, like 2D EMHD, classical and MHD turbulence, is isotropic when 
allowed to freely decay. In the presence of a background fiow, classical turbulence remains 
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FIG. 3: Plots of the Rb = 300 solution in real configuration space at t = 0.2. For full details of 
the slices through the datacube, see the text. 

isotropic. Small-scale structures are advected along by any large-scale flow, whether or 
not that has a uniform background contribution. This effect has been attributed to local 



coupling in phase space. Numerical simulations of MHD turbu 
strongly anisotropic in the presence of a background field ( 17, 



ence have found it to be 
m, Shebalin et al. 1983, 



Oughton et al. 1998]). This has been attributed to the excitation of Alfven waves which 
preferentially propagate parallel to the external magnetic field and hinder the cascade process 
perpendicular to the external field. 

In 2D EMHD turbulence, recent numerical studies employing hyperdiffusivity (0, Q, 
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FIG. 4: A plot of energy against time for Rb = 10, 30, 100 & 300, increasing from left to right. 

Dastgeer et al. 2000, Dastgeer and Zank 2003]) have revealed anisotropic behaviour with 
the cascade strongly inhibited parallel to the background field. This can only be the result 
of asymmetry in the nonlinear spectral transfer process relative to the external magnetic 
field. In the context of local energy coupling in Fourier space, mediation by whistler waves 
has been proposed as the only way this asymmetry could be achieved {[6:, Dastgeer et al. 
2000]), by the mechanism detailed by Q]. Most recently, the numerical study of 21^, WH09] 
has confirmed the anisotropy of 2D EMHD with normal diffusivity. The spectrum of 2D 
anisotropic EMHD turbulence has also been shown to exhibit a linear relationship with an 
external magnetic field ([^, Dastgeer and Zank 2003]). 

In order to understand how hyperdiffusivity has affected previous studies we have intro- 
duced a background field into the governing equations as discussed above and calculated 
solutions for Rb = 100 at a spatial resolution of 256^ points in real configuration space. We 
present our results in Figure O /^From left to right, we show 2D energy spectra for Rb = 100 
with C = 0, 1, 2, 4 & 8. In the isotropic case with no background field, i.e. C = 0, energy is 
evenly distributed between x and y, as indicated by circular contours. In the case of C = 1 
we find energy transfer to larger k has been suppressed in the x direction, parallel to the 
background field. 3D decaying EMHD turbulence has become anisotropic in the presence 
of a uniform background field with normal diffusivity. The effect becomes more pronounced 
for C = 2 and particularly strong for C = 4 and C = 8. It is at these high values of C 
that the cascade is almost turned off in the direction parallel to the field and even strongly 
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FIG. 5: 2D Fourier power spectra at t = 0.2 for 3D EMHD turbulence in the presence of a 
background field. We show three 2D Fourier power spectra slices of the Bx datacube, at kx = 1 
(top row), ky = 1 (middle row) and k^ = 1 (bottom row). Across the columns we show the power 
spectra for C = 0, 1, 2, 4 & 8. 

suppressed in the perpendicular directions. In 2D this suppression has been attributed to 
excitation of whistler waves, which act to weaken spectral transfer along the direction of 
propagation (j?!, Dastgeer and Zank 2003]). In Figure [6] we show a slice through the real 
configuration space field at 2; = 0.5. For C = 0, the field is isotropic, but as the value of 
C is increased, structures stretch in the x direction corresponding to increasingly inhibited 
energy transfer in the x direction but not in y. We find the same for x compared to z. At 
C = 8, the structures have remained almost the same as at t = since the cascade in all 
directions has been so strongly inhibited. 

jsj performed 3D simulations of EMHD turbulence. They introduced a uniform back- 
ground field of strength comparable to that of the fluctuating freely decaying field (i.e. 
C = 1). The authors employed a hyperdiffusivity of 3 and found scale-dependent anisotropy. 
Our simulation at C = 1 supports their result. Further, our simulations imply the linear 
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FIG. 6: Slices of real configuration space Bx fields at t = 0.2 for 3D EMHD turbulence in the 
presence of a background field. From left to right, we show the fields for C = 0, 1, 2, 4 & 8. In all 
cases z = 0.5. 

relationship between EMHD turbulence and strength of external magnetic field found by 7|] 
in 2D can be extended to 3D. 

Note that ([1]) is scale invariant, i.e. it is possible to apply the equation over the whole 
of a system, or just to a small section, with Rb unchanged. A very small box then will see 
the large-scale field as a background field, and therefore the smallest scales in the system, 
for example a neutron star, should be anisotropic. 

C. Intermediate-scale initial conditions 

In 2D classical turbulence, the exchange of energy and enstrophy Q is coupled in Fourier 
space according to 

^ - -e^ (6) 

dt~ " dt' 

hence energy injected at intermediate scales experiences a transfer to both higher and lower 
wavenumbers in order to satisfy this coupling and simultaneously conserve energy and en- 

Kraichnan 



strophy. This is the inverse cascade of energy to lower k (larger scales) ( [161 
1967]). In MHD turbulence, energy and magnetic helicity are coupled in the same way 
and an inverse cascade occurs in order to simultaneously conserve these two quadratic ideal 



invariants. In driven 2D EMHD tubulence, the inverse cascade found by 18| has been pro- 



posed to be the result of simultaneous conservation of energy and mean squared magnetic 



potential, or anastrophy. Our recent simulations ([2l|, WH09]) have shown an inverse cas- 
cade occurs in freely decaying 2D EMHD turbulence. Whilst it has a considerably different 
scaling nature to driven 2D EMHD turbulence, we believe it is also the result of simultane- 
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Fourier power spectra 
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FIG. 7: Time evolution of the Fourier power spectra for e = 2.5 -10 . Shown are power spectra 
at t = 0, 0.01, 0.02, 0.03 and 0.04. No further inverse cascase develops after t = 0.04. 

ous conservation of energy and anastrophy. In 3D then, the existence of an inverse cascade 
is immediately under question since anastrophy is no longer conserved. We inject energy 
over the wavenumber range 10 < /c < 20 as detailed above and evolve the magnetic field to 
assess what inverse cascade there is, if any. 

In this case, we have found that simulations with normal diffusivity are unable to pro- 
duce an inverse cascade. Since the values of Rb available to us are considerably less than 
Rb = 1000, where decaying 2D EMHD turbulence displays an inverse cascade, this is not a 
surprising result. For the runs in this section we therefore introduce a hyperdiffusivity, so 
that the large scales see a diffusivity equivalent to Rb > 1000, whereas the small scales are 
much more strongly damped, so that the computations can be done at all. As noted above, 
using a hyperdiffusivity will of course disrupt features like the precise shape of the spectrum 
at large k, including the presence or absence of a dissipative cutoff. However, since here we 
are interested in the behaviour at more modest k, in particular the possibility of transferring 
energy from intermediate to small wavenumbers, these distortions at large k probably do 
not have that much effect on the results. 

The precise form of hyperdiffusivity is simply to raise the power of the diffusion operator, 

^ = -V X [J X B] - e (WYB, (7) 

where e is the new diffusion coefficient, the equivalent of Rb^, but only as seen by the largest 
scales. 
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Figure [7] shows spectra at t = 0, 0.01, 0.02, 0.03 and 0.04. The spectra show that energy 
is transferred to < 10 in a very weak inverse cascade. The spectral peak is slowly shifting 
to < 10 but not maintaining the same amplitude. Some energy has also been transferred 
to A; > 20. In 3D then, any inverse cascade is considerably weaker than in 2D, presumably 
since anastrophy is no longer conserved. 



IV. CONCLUSIONS 



We have investigated the nature of decaying 3D EMHD turbulence with normal diffusivity 
and compared it with classical and MHD turbulence and studies of 2D and 3D EMHD with 
hyperdiffusivity. We have found 3D EMHD turbulence experiences an isotropic forward 
cascade of energy to higher wavenumber (smaller spatial scales) asymptotically approaching 
Ek oc A;~^ with increasing Rb (inversely proportional to a dissipation coefficient). This is 
in agreement with the original theoretical prediction of lOj, suggesting freely decaying 3D 
EMHD turbulence is mediated by whistler waves. Unlike the theoretical prediction, we have 
found there is no dissipative cutoff at the predicted wavenumber k Rb and argue this 
is consistent with non-local coupling in Fourier space, the most important result of this 
paper and consistent with our recent results for 2D EMHD turbulence. Hyperdiffusivity has 
previously clouded this issue and introduced an artificial cutoff and the unwanted bottleneck 
effect. We have also found that fully developed 3D EMHD turbulence appears to be strongly 
structured, retaining a similarity to the initial field at late time, very much unlike classical 
or MHD turbulence. 

3D EMHD turbulence with normal diffusivity has been found to display scale-dependent 
anisotropy in the presence of a uniform background field, in agreement with previous studies 
employing hyperdiffusivity. Our results support previous studies which found the strength 
of the anisotropy is linearly related to the external field strength. Further, strong fields 
effectively halt the cascade in the parallel direction and considerably inhibit the cascade in 
all other directions. 

Finally, we have discovered that decaying EMHD turbulence yields a weak inverse cas- 
cade, at least partially transferring magnetic energy from intermediate to large lengthscales. 
The possibility of inverse cascades may have implications for the magnetic fields of neutron 
stars, where the proto-neutron star that emerges from a supernova explosion may well have a 
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primarily small-scale, disordered field. A Hall-induced inverse cascade may then be a mech- 
anism whereby it acquires a large-scale, ordered field. We note though that the electron 
number density n is strongly depth-dependent in neutron stars, which turns out to interact 
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Vainshtein et al. 2000, Hollerbach 



with the Hall effect in a highly nontrivial way 
and Riidiger 2004]). It is likely therefore that both forward and inverse cascades will be 
rather different in real neutron stars than they are in isotropic, Cartesian box models such 
as here. Future work will consider EMHD in more realistic, stratified spherical-shell models. 
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